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SUMMARY 


[  Thq'present  experimental  ("study  is  concerned  with  the  effect  of  high- 
velocity  fluid  flow  on  the  static  and  dynamic  characteristics  of  a  simply  sup¬ 
ported  pipe.  The_ frequency  variation  of  the  bending  vibrations  as  a  function 
of  fluid  velocity^  as  predicted  by  G.  W.  Housner,  (was  verified  experimentally. 
In  addition,  the  critical  fluid  velocity  at  which  'the  system  becomes  statically 
unstable  was  also  verified.\  \ 

H.  Nf'  J 


INTRODUCTION 


High-performance  liquid -propellant  launch  vehicles  require  the  rapid 
transfer  of  large  quantities  of  fluid  from  tanks  to  pumps  and  engines  through 
relatively  thin-walled  pipes.  A  study  of  the  influence  of  the  resulting  high- 
velocity  fluid  on  the  static  and  dynamic  characteristics  of  the  pipes  in  launch 
vehicles  is  therefore  necessary. 

Initial  work  in  this  field  was  motivated  by  vibration  problems  encountered 
in  the  transport  of  crude  oil  in  pipes  as  well  as  possible  problems  in  the  feed 
lines  of  large  water  turbines.  Ashley  and  Haviland  (ref.  l)  were  the  first  to 
investigate  the  bending  vibrations  of  a  simply  supported  pipe  containing  a 
flowing  fluid.  Housner  (ref.  2)  subsequently  derived  the  equation  of  motion 
for  the  same  case  more  completely  and  developed  an  equation  relating  the  fun¬ 
damental  bending  frequency  of  a  simply  supported  pipe  to  the  velocity  of  the 
internal  flow  of  the  fluid.  Housner  also  showed  that,  at  a  certain  critical 
velocity  (similar  to  the  Euler  buckling  formula  for  columns),  a  statically 
unstable  condition  could  exist.  Long  (ref.  5)  used  an  alternate  solution  to 
Housner1 s  equation  of  motion  for  the  simple  end  conditions  and  also  treated  the 
fixed-free  end  conditions.  In  addition,  Long  attempted  to  confirm  the  analysis 
by  comparison  with  experiment.  The  experimental  results  were  rather  inconclu¬ 
sive  since  the  maximum  fluid  velocity  available  for  the  test  was  low  and  the 
change  in  the  bending  frequency  was  very  small.  The  maximum  velocity  obtained 
by  Long  was  only  20  percent  of  the  critical  velocity  which  resulted  in  a 


maximum  change  in  the  tending  frequency  of  only  3-2  percent.  Other  treatments 
of  the  subject  were  made  by  Ta  Li  (ref.  4),  Niordson  (ref.  5) ,  and  Benjamin 
(ref.  6) . 

It  is  the  purpose  of  this  paper  to  report  the  results  of  an  experiment  to 
(a)  study  the  effect  of  fluid  velocity  on  the  fundamental  bending  frequency  of 
a  simply  supported  pipe  and  (b)  determine  the  existence  of  a  statically 
unstable  condition  as  predicted  in  reference  2. 

A  number  of  engineering  students,  while  participating  consecutively  in 
the  cooperative  work  study  program  at  Langley  Research  Center,  conducted  the 
test  program  under  the  general  supervision  of  the  second  author.  The  students 
who  planned  the  experiment,  selected  and  assembled  the  test  apparatus,  and 
conducted  the  tests  were  Messrs.  Harold  L.  Dodds,  Jr.,  Richard  A.  Crocker, 
Charles  H.  Pox,  Jr.,  James  E.  Gardner,  and  A.  David  Luckey. 
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SYMBOLS 


amplitude  coefficient  in  series  expression  (eq.  (2))  for 

k  =  1,  2,  3,  •  •  • 

amplitude  coefficient  in  series  expression  (eq.  (2))  for 
n  =  1,  2,  3,  .  .  . 

modulus  of  elasticity,  lb/ft^ 

moment  of  inertia,  ft^ 

integer 

length,  ft 

mass  of  pipe  plus  fluid  per  unit  length,  slugs/ft 

integer 

time,  sec 

average  velocity  of  the  fluid  relative  to  x-direction,  ft/sec 
coordinate  measured  along  pipe 
coordinate  measured  perpendicular  to  pipe 
mass  of  fluid  per  unit  length,  slugs/ft 

angular  velocity  of  pipe  relative  to  its  length,  rad/sec 
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Subscripts: 


c  critical 

i  mode  of  vibration,  i  =  1,  2,  3,  .  .  . 

ANALYSIS 


The  equation  of  motion  of  a  uniform  pipe  containing  flowing  fluid  is 


El 


^x4 


_PV2  ^5.  _  2pV  ■■ 

ax2  St  Sx 


-  m 


St2 
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which  was  derived  in  reference  2  by  using  elementary  beam  theory.  The  coordi¬ 
nate  system  is  shown  in  figure  1.  The  three  terms  on  the  right  side  of 


Figure  1.-  Schematic  diagram  of  experimental  setup. 


equation  (l)  represent  the  lateral  inertial  forces  which  act  on  the  pipe.  The 
first  term  is  the  force  caused  by  the  change  in  direction  of  velocity  of  the 
fluid  due  to  the  curvature  of  the  pipe.  The  second  term  is  the  force  associ¬ 
ated  with  the  Coriolis  acceleration  which  occurs  because  the  pipe  has  an 
angular  velocity  with  respect  to  its  length  while  the  fluid  has  a  linear  veloc¬ 
ity  with  respect  to  the  pipe.  The  last  term  represents  the  force  which  is 
related  to  the  vertical  acceleration  of  the  pipe.  For  simplification,  struc¬ 
tural  damping  has  been  omitted  in  equation  (l). 

For  the  case  of  free  vibrations,  the  general  solution  of  equation  (l)  for 
a  simply  supported  pipe  with  normal  modes  z p  is  given  in  reference  2  as 
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In  order  to  evaluate  the  coefficients  in  this  equation,  equation  (2)  is  sub¬ 
stituted  into  equation  (l)  and  the  resulting  terms  are  grouped  according  to 
their  trigonometric  constituents.  Equating  the  coefficients  of  like  trigono¬ 
metric  terms  to  zero  in  order  to  satisfy  equation  (l),  the  following  set  of 
equations  for  n  =  1,  2,  3>  .  .  .  and  k  =  1,  2,  3,  .  .  .  are  obtained: 
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Supressing  all  the  modes  except  the  first  two,  the  simultaneous  solution  of 
equations  (3)  and  (4)  results  in  the  following  frequency  equation  for  steady- 
state  vibrations: 


32pVoai 
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Equation  (5)  gives  the  relation  between  V  and  <pj_  for  the  pipe. 
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By  letting  toy  -»  0  for  the  case  i  =  1,  equation  (5)  yields  a  critical 
velocity  at  which  a  static  divergence  occurs. 


Vc  = 


Elm2 


(6) 


As  pointed  out  in  reference  2,  replacement  of  pV2  in  equation  (6)  hy  P 
provides  an  equation  identical  to  the  Euler  buckling  formula  for  columns  where 
P  represents  the  axial  load. 


EXPERIMENTAL  INVESTIGATION 


The  apparatus  and  testing  procedure  which  were  used  to  study  the  effect  of 
high-velocity  fluid  flow  on  the  bending  vibrations  and  stability  of  a  simply 
supported  pipe  are  described  in  the  following  section. 


Apparatus 

A  schematic  representation  of  the  experimental  setup  is  shown  in  figure  1. 
The  test  pipe  was  constructed  of  aluminum  alloy  (6o6l~T6)  with  a  nominal  out¬ 
side  diameter  of  1.00  inch,  a  wall  thickness  of  0.065  inch,  and  a  length  of 
10  feet.  The  pipe  was  mounted  horizontally  with  the  downstream  end  supported 
by  two  ball  bearings  which  allowed  unrestrained  horizontal  translation  and 
rotation  of  the  downstream  end.  The  upstream  end  was  connected  to  a  swivel 
joint.  This  type  of  support  system  provided  an  accurate  simulation  of  a  beam 
with  pinned-pinned  boundary  conditions.  The  ball  bearings  at  the  end  of  the 
pipe  prevented  axial  stresses  during  lateral  deformations.  A  strain  gage  was 
mounted  in  the  center  of  the  pipe  at  the  location  of  the  maximum  bending  moment 
from  which  the  frequency  and  damping  of  the  bending  vibrations  could  be 
determined. 

The  average  velocity  of  the  fluid  (water)  was  measured  by  means  of  an 
orifice  plate  located  approximately  4  feet  upstream  from  the  swivel  joint. 

The  instrumentation  for  measuring  the  fluid  velocity  consisted  of  the  orifice 
plate,  a  differential-pressure  transducer,  and  a  galvanometer  recorder.  This 
equipment  was  accurate  to  within  2  percent  (rms  value)  of  full-scale  response. 
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The  fluid  for  this  experiment  was  obtained  from  a  high-pressure  water 
source  located  at  the  Langley  landing  loads  track.  This  facility,  capable  of 

delivering  750  gallons  per  minute  at  a  maximum  pressure  of  3200  lb/in2,  is 
described  in  detail  in  reference  7- 


Testing  Procedure 

In  order  to  initiate  the  free  vibrations  of  the  pipe  at  various  fluid 
velocities,  a  method  of  imparting  a  small  disturbance  to  the  pipe  was  necessary. 
This  disturbance  was  accomplished  by  attaching  a  cord  at  the  midpoint  of  the 
pipe  which  could  be  manually  operated  during  the  test.  The  procedure  for  con¬ 
ducting  a  test  was  to  disturb  the  pipe  by  plucking  the  cord  vertically,  to 
record  the  response  of  the  pipe,  and  to  increase  the  fluid  velocity  continu¬ 
ously  until  instability  of  the  pipe  was  achieved. 


DISCUSSION  OF  RESULTS 


This  section  includes  both  the  calculated  and  experimental  results  and 
also  some  remarks  concerning  the  "damping"  characteristics  of  the  system. 

The  values  of  the  constants  used  in  calculating  the  value  of  Vc  from 
equation  (6)  and  the  relation  between  o>]_  and  V  from  equation  (5)  are  as 
follows: 


p  =  8.00  X  10 "3  slugs/ft 

E  =  1.44  X  109  lbs/ft2 
I  =  1 . 04l6  X  10"6  ft^ 
m  =  1.49  X  10-2  slugs/ft 
Z  =  10.5  ft 

A  departure  from  an  ideally  uniform  beam  was  required  for  the  upstream 
end  of  the  pipe  in  order  to  accommodate  part  of  the  swivel  joint  (fig.  2). 

The  joint  provided  about  0.5  foot  of  rigid  heavy  pipe  for  an  overall  test 
length  of  10.5  feet.  The  effect  of  this  short  stiff  section  (compared  with 
the  thin-walled  aluminum  tube)  on  the  critical  velocity  is  considered  negligi¬ 
ble  since  the  curvature  of  the  pipe  is  a  maximum  at  its  center  and  must 
approach  zero  at  the  ends  in  order  to  satisfy  the  boundary  conditions.  On  the 
other  hand,  the  effect  of  the  swivel  joint  on  the  frequency  of  the  pipe  could 
be  more  pronounced  because  of  the  concentrated  mass  at  the  end  of  the  pipe. 
However,  comparison  of  the  calculated  frequency  for  a  uniform  pipe  10. 5  feet 
long  with  the  experimental  frequency  for  the  test  pipe,  which  includes  the 
deviation  resulting  from  the  swivel  joint,  indicates  no  appreciable  effect  on 
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Figure  2.-  Swivel  joint  at  upstream  end  of  test  pipe. 
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the  natural  frequency.  Consequently,  the  fundamental  equations  in  refer¬ 
ence  2  were  not  rederived  to  include  the  possible  effects  caused  by  the 
swivel  joints. 


The  calculated  results  from  equations  (5)  and  (6)  for  the  fundamental 
mode  of  vibration  and  the  experimental  data  for  two  identical  pipes  are  shown 
in  figure  3*  The  specific  experimental  values  pertaining  to  figure  3>  which 
have  been  normalized  by  the  theoretical  values  of  cd  at  V  =  0  and  Vc 
^oy-o  =  28.4  rad/sec  and  Vc  =  129-5  ft/sec  appear  in  the  following  table. 

The  overall  damping  of  the  system  at  each  data  point,  as  determined  from  the 
strain-gage  records,  also  appears  in  the  table: 


Pipe 

v/vc 

s 

<? 

II 

O 

Log  decrement 

*0 

1.042 

0.179 

]_ 

•332 

•919 

.308 

•  595 

.849 

.409 

•753 

.662 

.967 

0 

Not  obtained 

Not  obtained 

.167 

1.053 

0.329 

•  354 

.958 

•  337 

2 

.543 

.919 

.629 

.689 

.676 

.865 

•  752 

•739 

1.098 

.985 

0 

Not  obtained 

*For  pipe  1  at  v/vc  =  0,  the  data  were 

obtained  when  the  pipe  was  empty;  however, 
the  value,  =  1.042,  has  been  cor¬ 

rected  to  include  the  mass  of  the  fluid  for 
a  full  pipe . 
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Figure  3.-  Comparison  of  theoretical  and  experimental  results. 


As  can  be  observed  in  figure  3>  a  good  correlation  between  theory  and 
experiment  is  achieved.  The  maximum  experimental  error  relative  to  the  the¬ 
oretical  curve  is  approximately  7*5  percent. 

As  predicted  by  the  theory,  the  system  became  unstable  which  resulted  in 
a  permanent  deformation  of  the  pipe.  Divergence  was  achieved  with  both  pipes; 
however,  because  of  an  instrumentation  failure  during  the  first  test,  only  the 
divergence  velocity  for  the  second  pipe  was  recorded  which  was  127-5  ft/sec. 
The  experimental  error  of  the  divergence  velocity  relative  to  the  predicted 
critical  velocity  is  1.5  percent.  Pictures  of  the  divergence  phenomenon  are 
shown  in  sequence  in  figure  4.  Figure  4(a)  shows  the  curvature  of  the  pipe  at 
V  =  120  ft/sec  which  is  just  before  the  pipe  diverged.  Figure  4(b)  shows  the 
curvature  of  the  pipe  at  the  divergence  velocity  of  127-5  ft/sec.  Figure  4(c) 
shows  the  permanent  deformation  of  the  pipe  after  the  test  was  completed. 

The  log  decrement,  as  shown  in  the  table,  is  a  measure  of  the  total 
damping  of  the  system  at  each  data  point .  These  values  were  obtained  by 
taking  the  natural  logarithm  of  the  ratio  of  successive  amplitudes  for  the 
first  completed  cycle  of  each  strain-gage  record.  The  log  decrement  value 
when  the  pipe  is  empty  can  be  considered  as  a  measure  of  the  structural 
damping  of  the  pipe  and  its  end  supports  (ref.  3).  This  value  was  0.179  for 
the  first  pipe. 

As  may  be  observed  in  the  table,  the  natural  frequency  of  the  pipe 
decreases  as  the  fluid  velocity  increases  while  the  log  decrement  values 
increase.  It  should  be  noted  that  the  log  decrement  increase  is  due  to  the 
combined  effect  of  the  frequency  and  fluid  velocity,  although  the  values  dif¬ 
fer  in  some  cases  for  the  two  pipes. 
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CONCLUDING  REMARKS 


The  present  investigation  has  been  concerned  with  an  experimental  study 
of  the  effect  of  fluid  velocity  on  the  natural  frequency  of  the  free  vibra¬ 
tions  of  a  simply  supported  pipe.  Tests  were  conducted  with  two  identical 
pipes  both  of  which  verified  the  predicted  variation  of  the  fundamental 
bending  frequency  with  fluid  velocity.  The  predicted  critical  fluid  velocity 
for  which  the  system  becomes  statically  unstable  was  also  verified  experimen¬ 
tally.  From  the  results  presented,  it  appears  that  the  theory  of  G.  ¥.  Housner 
is  adequate  for  predicting  the  effect  of  fluid  velocity  on  the  natural  fre¬ 
quency  of  simply  supported  pipes  as  well  as  the  particular  fluid  velocity  at 
which  a  static  instability  exists. 


Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Langley  Station,  Hampton,  Va.,  March  15,  1965- 
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